BCS-to-BEC crossover from the exact BCS solution 
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The BCS-to-BEC crossover, as well as the nature of Cooper pairs, in superconducting and Fermi superfluid 
media is studied from the exact ground state wavefunction of the reduced BCS Hamiltonian. As the strength 
of the interaction increases, the ground state continuously evolves from a mixed-system of quasifree fermions 
and pair resonances (BCS), to pair resonances and quasibound molecules (pseudogap), and finally to a system 
of quasibound molecules (BEC). A single unified scenario arises where the Cooper-pair wavefunction has a 
unique functional form. Several exact analytic expressions such as the binding energy and condensate fraction 
are derived. We compare our results with recent experiments in ultracold atomic Fermi gases. 
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The nature of Cooper pairs in the BCS-BEC crossover has 
regained attention due to the observation of a large fraction 
of preformed fermion pairs on the BCS side of the Feshbach 
resonance in ultracold atomic Fermi gases Q|. While several 
theoretical explanations were proposed Q|, the interpretations 
are still controversial. The root of the controversy is under- 
standing what represents a Cooper pair in a correlated Fermi 
system, a concept not clearly defined in the literature. This 
paper discusses this concept at the very foundational level and 
proposes a qualitatively different scenario of the BCS-BEC 
crossover, based on the exact solution to the BCS Hamiltonian 
01- Only in the extreme BEC limit does this new scenario and 
Leggett's naive ansatz 0] become identical. 

While the superconducting and Fermi superfluid thermody- 
namic states represent a mixed-system of quasifree and pair- 
correlated fermions, the molecular BEC which arises in the 
dilute and strong coupling limit has all fermions bound into 
pairs forming a unique macroscopic quantum state. It is by 
now well accepted in which sense these states represent a 
Bose-Einstein (BE) condensation. What defines a BE con- 
densation in an interacting TV-particle system is spontaneous 
gauge symmetry breaking, or phase coherence, of its quan- 
tum state (a concept that strictly applies in the thermodynamic 
limit (TL)). Yang |5] provided a detailed mathematical char- 
acterization of this phenomenon which manifests itself as off- 
diagonal long-range order (ODLRO) or, equivalently, by the 
existence of an eigenvalue of order N in a reduced density 
matrix p n , where n is the number of particles forming the 
smallest unit that condenses. 

The above definition of a BEC does not imply the naive 
view of a many-body state of quantum objects with identical 
internal wavefunctions. Indeed, we will argue that a current- 
carrying mean-field ground state (GS) of a uniform supercon- 
ducting or Fermi superfluid iV-particle system is of the form 

• ■ • ,x N ) = A [(pi(xi,x 2 ) ■ ■ ■ <f>N/2(xtf-l,Xtf)] , (1) 
with Xj = (tj ,<Jj), antisymmetrizer A, and the pair state 

- (rj+l-j) 

4> a (xi,Xj) = e q ' ' tp a (n - rj) x{o-i, o-j) , (2) 

where \ is a spin function of a certain symmetry, q is the 
pair center-of-mass momentum, and <p a {r) the internal wave- 



function which may represent either a quasimolecular reso- 
nant state or a scattering state (i.e., a mixed state), depending 
upon the strength of the interaction between particles. 

For pedagogical reasons we will recall the one-Cooper-pair 
problem and then address the question of what happens when 
we add more pairs (q = in the GS). The Cooper-pair solu- 
tion can be recovered from the wavefunction ([Q by assuming 
that N — 2 fermions are free, filling the lowest momentum 
states k (of energy e k ) up to the Fermi momentum kp, while 
only an additional fermion pair (with momenta k > kp) feels 
an attractive s-wave interaction in the spin singlet channel 

i*>c = e ^rtAi w ■ (3) 

The role of the Fermi sea, \F ) = IIfc<fc F c kT c -kj0)' ls t0 
Pauli -block states below the Fermi energy ep. Assuming that 
the attractive pairing interaction G < is constant around the 
Fermi energy, the eigenvalue E is always negative implying 
that the Cooper pair is bound regardless of the strength of the 
attractive interaction. The Fermi sea is then unstable against 
the formation of bound pairs of electrons. 

What happens when the pairing interaction also affects 
electrons in the Fermi sea? The answer to this question is 
the BCS theory whose canonical form considers a state of the 
form Q with all identical internal wavefunctions 

P M |BCS) = (At) M |0) , At = E -4 TC L ki , (4) 

k 

where Pm is the projector onto the state with M pairs, and i?k, 

Mk are the BCS amplitudes (J = \ (\ T -^-Sgi—) 

with A the gap and fj, the chemical potential. Since the struc- 
ture of the BCS pair is averaged over the whole system, it 
does not characterize a Cooper pair in the superconducting 
medium except in the extreme strong-coupling and dilute lim- 
its where all pairs are bounded and non-overlapping. Some- 
times, the structure of the Cooper pair is associated with the 
pair-correlation function (BCS|cj c .j.c_ k jBCS) = u k v k de- 
scribing the pair correlation among fermions of the same pair 
as well as the exchange between fermions of different pairs. 



What is the nature of a Cooper pair in a Fermi superfluid 
or superconducting medium? To address this question we will 
use the exact solution to the reduced BCS Hamiltonian 
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proposed by Richardson 40 years ago 0>laU- H involves all 
terms with time-reversed pairs (k |, — k J.) from a contact 
interaction. It is consistent with an effective single-channel 
description of the BCS-BEC crossover theory 0] in terms of 
a zero-range potential. Realistic finite-range interactions pro- 
duce qualitatively similar results along the crossover Q| . 

For simplicity we will consider TV = iVj + N± spin- 1/2 
(i.e., 2-flavor) fermions in a three-dimensional box of volume 
V with periodic boundary conditions, interacting through an 
attractive constant (s-wave-singlet-pairing) potential. (Other 
pairing symmetries can also be accommodated.) Exact N = 
2M + v particle eigenstates of H can be written as 
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where \u) = \v\, v% ■ ■ ■ , vl) is a state of v unpaired fermions 
(v = Sk^k, with v\t = 1,0) defined by c_ k ^c k1 - \v) = 0, 
and ?ik \ v) = vv. \ v )- L is the total number of single particle 
states. The GS |* ) is in the v = (7V T = N{) sector. 

Each eigenstate \*f?) is completely defined by a set of M 
(in general, complex) spectral parameters (pair energies) Et 
which are a solution of the Richardson's equations 
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and the eigenvalues of H are £ = ^ k £k^k + YlfiLi Ee _ 
A crucial observation is that if a complex Eg satisfies 0, its 
complex-conjugate E% is also a solution. Thus, restores 
time-reversal invariance. The ansatz (jfji is a natural general- 
ization of the Cooper-pair problem without an inert Fermi sea, 
and with all pairs subjected to the pairing interaction. The pair 
structure in l|6} is similar to the Cooper pair in 0, and the 
many-body state has the form Q with all pairs different. 

Since we are concerned with uniform bulk Fermi systems, 
we are interested in the TL (i.e., N, V — > oo with p = N/V = 
kp/3TT 2 = constant). This limit was studied by Gaudin 0] 
when the energy spectrum G [— u>, u>) is bounded, and v = 
0. Eqs. reduce to the well-known gap equation 
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(e - pf + A 2 
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where 17(e) represents the density of states. In our case 
fi € [0,oo] is unbounded with g(e) = m 3/2 e 1/2 /(V2tt 2 H 3 ) 
for £k = H 2 k 2 /2m. Due to the absence of an upper cutoff, 
Eq. is singular and the TL becomes a subtle mathematical 
procedure whose solution will be presented here 1 10]. 
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TABLE I: Analytic expressions for selected values of 77 = l/fcj?a s ; 
x is the root of Pi , i.e. Pi (x) = 0, and E(y) is the complete elliptic 

integral of the second kind. Note that 5ir 2 £o(r] = 0)/k% = p(r) = 
0) « 0.590606. 



There are two ways to regularize the problem: One can ei- 
ther introduce an energy cutoff in the integrals or one can can- 
cel the singularities with physical quantities whose bare coun- 
terpart diverges in the same way Q. For this problem, Eq. 
(|5}, the physical quantity is the scattering length a s given by 
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The non-singular gap equation (after integration fiHl ) is 



kpa s 



A 2 Pi 



vV + A 2 



, (10) 



where energies are now in units of ep = h 2 k F /2m and 
lengths in units of £p = 1/kp. Pp(x) is the Legendre func- 
tion of the first kind of degree j3. The equations for the con- 
servation of the number of pairs M 
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and GS energy density, for arbitrary strength 77, 
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do not need regularization, these are exact results. Indeed, 
for a given 77, the chemical potential p and gap A need to be 
determined self-consistently from Eqs. JlOi and M 11 (Fig.0. 
Then, the GS energy can be computed as a function of density 
using Eq. Mil . It shows no phase segregation. The exact 
binding (or condensation) energy per electron pair £g is (see 
Tabled 
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The complete solution of Eqs. @ in the TL 0] amounts to 
determining (for a given 77) the set of pair energies Eg which, 
for the GS, form a single open, continuous, and reflection- 
symmetric arc r with extreme points Ep = 2(p ± iA). Here, 
we simply present the results without derivation. The equation 
for r is 
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FIG. 1: Ground state crossover diagram displaying the different 
regimes as a function of rj. A few arcs T, whose extremes are 
E = 2(/i ± iA), are displayed. The dashed line corresponds to 
the Feshbach resonance rj — 0. Open circles represent the values of 
E for which Cooper-pair wavefunctions are plotted in Fig. [5] Insets 
show the behavior of the gap A and chemical potential ji. 
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where z= y {E — p) + A 2 . 

Figure^shows the BCS-BEC crossover diagram and sev- 
eral arcs corresponding to particular values of r). Since a 
crossover diagram is not a phase diagram, i.e. there is no 
symmetry order parameter or non-analyticities sharply differ- 
entiating the regions, boundaries are in principle arbitrary. 

Here we have adopted the following criteria: The geome- 
try of the arcs T serves us to establish a criterium for defining 
boundary regions in the crossover diagram. In the extreme 
weak-coupling limit, rj — > — oo, the pair energies are twice 
the single particle energies (Et — > 2et). As soon as the in- 
teraction is switched on, i.e. kpa s is an infinitesimal negative 
number, a fraction f of the pair energies close to 2ji become 
complex, forming an arc T in the complex plane. The fraction 
1 — f of fermion pairs below the crossing energy 2e c of T with 
the real axis have real pair energies. They are still in a sea 
of uncorrelated pairs with an effective Fermi energy e c . The 
dark grey region labelled BCS, which extends from rj = — oo 
to rj = 0.37, is characterized by a mixture of complex pair en- 
ergies with a positive real part and real pair energies. rj = 0.37 
is the value at which all pair energies are complex, i.e. f = 1, 
and the effective Fermi sea has disappeared. Within the BCS 
region we plotted the arcs for rj = — 1 (solid line), with a 



fraction f = 0.35, and rj = (dashed line) with f = 0.87. 
The pseudogap region P, indicated in light grey, extends from 
rj = 0.37 to rj = 0.55 where p, = 0. Within this region the 
real part of the pair energies changes from positive to neg- 
ative, and P describes a mixture of Cooper resonances and 
quasibound molecules. The BEC (white) region, rj > 0.55, is 
characterized by all pair energies having negative real parts, 
i.e. all pairs are quasibound molecules. As rj increases fur- 
ther, r tends to an almost vertical line with Kc(E) ~ 2/x, and 
-2A < Im(E) < 2A. 

If the literature is not clear about the size £ of a Cooper 
pair, it is partly because it is not clear what a Cooper pair 
is. Pippard, in his nonlocal electrodynamics of the super- 
conducting state, introduced the characteristic length £o by 
using an uncertainty-principle argument involving only elec- 
trons within a shell of width 2A around the Fermi energy. 
The resulting quantity, named Pippard's coherence length, 
£o = 2/ (ttA), is usually associated to £. On the other hand, an 
acceptable definition could be £ = \J ((p\r 2 \tp) , (<p\<p) = 1. 



From Eq. 0, the Cooper-pair wavefunction is 
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with A 2 = Im( v /s72)/2 7 r£3 , ipg = C/(2e k - E), and C 



being a normalization constant. Thus, £g = l/Im(\/£')- In 
the weak-coupling BCS limit (A <C ji ~ 1), when E = Ef, 
we get £e = ^S,o/V2. On the other hand, in the same limit, if 
one uses = Cp^k^k, one gets £p = £e/2, and if one uses 
^k CS = CecsWitk, one gets £ B cs = V 21 / 2 - 
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FIG. 2: Different definitions of Cooper-pair sizes £, and condensate 
fraction A and f as a function of the interaction strength rj. The hori- 
zontal solid line represents the interparticle distance r s — \J 97r/4. 

Analytic expressions (using dimensional regularization) for 
£p and £bcs can also be obtained for arbitrary coupling 
strength. Figure [2] displays the different sizes, all labelled 
by £, as a function of rj. Clearly, in there is no unique 
pair size £e but a distribution of sizes; Fig. |2]shows in solid 
line the smallest size in T corresponding to Ep. On the other 
hand, there is a unique pair size for ipp and (^bcs; the dashed 
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and dotted lines correspond to the £p and £bcs sizes respec- 
tively. Notice the relation between sizes and the interparticle 
distance r s = v/9?r/4 which is indicated as a full horizon- 
tal line. While £g and £p increase for negative r\ values, and 
eventually diverge when r\ — * — oo, £bcs tends to a constant 
value, related to the interparticle distance, showing its essen- 
tially uncorrelated nature. The fact that £p < £e/2 in the 
overlapping BCS region is a clear demonstration that £p mea- 
sures the mean distance between time-reversed pairs irrespec- 
tive of the Cooper pair they belong to. For non-overlapping 
pairs (BEC region) both sizes converge to the same values. 
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FIG. 3: Cooper-pair sizes £e along the arcs T for the values of r\ 
depicted in Fig. 1. 

Following Yang H, ODLRO in p 2 (n t r 2 J. | r 3 1 r 4 J.) -► 
i/)*(rit, r 2 |)</>(r 3 T, r 4 |) may be used to define 
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as a measurement of the condensate fraction. Figure |2] shows 
A and the fraction f of (Cooper) pairs in the arc, that is, the 
fraction of complex spectral parameters. Although A can qual- 
itatively describe the fraction of correlated pairs, the fraction f 
gives the value 1 at the BCS-pseudogap boundary (?; = 0.37) 
while A = 1 for r\ — > oo. We note that at resonance (rj = 0), 
f = 87% of the fermions form Cooper pairs (A « 0.7). These 
results provide a qualitative explanation of the experiments in 
(jl without resorting to a projection method | 2|. Close to res- 
onance on the BCS side, a fraction f ~ 80% of the atoms form 
Copper pairs with sizes comparable to r s . Those atom pairs 
are efficiently transformed into quasimolecules during a rapid 
magnetic field ramping across the resonance giving rise to the 
molecular condensate fractions observed experimentally. 

One may argue that selecting the smallest pair from each arc 
to compare (in Fig0 Cooper-pair sizes along the crossover 
might not be representative of the Cooper-pair distribution 
within each arc. In Fig. [3] we show the internal variations 
of the pair sizes within the arcs. Although in the BCS re- 
gion there is a distribution of sizes from the smallest pair with 
E = Ep in the extremes of the arc to an almost infinite size 



pairs when they are close to the real axis, already at resonance 
(rj = 0) most of the pairs have the same structure. This latter 
feature becomes more pronounced in the BEC region. 

The various definitions of Cooper-pair wavefunctions are 
depicted in Fig. [4] which compares <^BCs( r )' <Pp(t), an d 
<Pe(t) for interaction strengths which correspond to the BCS, 
Feshbach resonance, and BEC regions of Fig. \l\ Notice 
that while tpBCS* an d ¥p aie highly oscillatory in the weak- 
coupling region, this is not the case with ipp which always 
has an exponential form. Clearly, a single and unified pic- 
ture emerges for the crossover when using a many-body state 
such as Q: tpE evolves smoothly through the crossover as it 
should. It is important to mention that the three wavefunctions 
are exactly the same in the strong coupling limit rj — > +oo. In 
this limit the naive ansatz of Leggett | 4] and the GS Q coin- 
cide, becoming a simple Pfaffian state. 




FIG. 4: Cooper-pair wavefunctions for different r\. The upper and 
middle (Feshbach resonance) panels correspond to the BCS region, 
while the bottom one is in the BEC region. Except for the Cooper 
ifE case, the other two wavefunctions always vanish at r = 0. It is 
only in the limit r\ — > +oo that the three states exactly coincide. 

In summary, we studied the BCS-BEC crossover problem, 
as well as the nature of Cooper pairs in a correlated Fermi sys- 
tem, from the exact GS l^o) of the reduced BCS Hamiltonian. 
We have analytically determined its exact TL for the quadratic 
single-particle dispersion, and calculated several quantities of 
physical interest, such as the binding (or condensation) en- 
ergy and the condensate fraction. The validity of the present 
description is not restricted to integrable pairing Hamiltoni- 
ans, but rather the ansatz l^o). which is a natural extension 
of the Cooper problem, could be used as a variational state 
for more realistic interactions. The Cooper-pair wavefunc- 
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tion ipE has a fixed functional s-wave form that accommo- 
dates pair-correlated fermions close to the Fermi energy in 
the weak coupling limit as well as quasibound molecules in 
the BEC limit. It is free from the characteristic oscillations 
displayed by </?bcs and ipp related to exchange among pairs. 
The BCS region in the crossover diagram represents a mixture 
of Cooper pairs and quasifree fermions. It may be argued that 
the single-channel model is insufficient to describe the system 
close to resonance where the degrees of freedom associated 
to the molecular closed channel should be explicitly incorpo- 



rated. A derivation analogous to the one presented here can 
be pursued by using a recently proposed atom-molecule inte- 
grable model fl2ll which captures the essential features of the 
two-channel model foil . The structure of the Cooper pair in 
this new model is the same as in (j5Jl. 
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